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Abstract 

The rat sarcoma-extracellular signal regulated kinase mitogen-actlvated protein kinases pathway, one of the most ancient 
signaling pathways, Is crucial for the defense against Bombyx mori nucleopolyhedrovirus (BmNPV) infection. Sprouty (Spry) 
proteins can inhibit the activity of this pathway by receptor tyrosine kinases. We cloned and identified a new B. mori gene 
with a Spry domain similar to the Spry proteins of other organisms, such as fruitfly, mouse, human, chicken, Xenopus and 
zebrafish, and named It BmSpry. The gene expression analysis showed that BmSpry was transcribed in all of the examined 
tissues and In all developmental stages from embryo to adult. BmSpry also induced expression of BmNPV In the cells. Our 
results Indicated: (1) the knock-down of BmSpry led to increased BmNPV replication and silkworm larvae mortality; (2) over- 
expression of BmSpry led to reduced BmNPV replication; and (3) BmSpry regulated the activation of ERK and inhibited 
BmNPV replication. These results showed that BmSpry plays a crucial role In the antiviral defense of the silkworm both in 
vitro and in vivo. 
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introduction 

Sprouty (Spry) is a general inhibitor of receptor tyrosine kinases 
(RTKs), first identified in fibroblast growth factor (FGF)-stimulat- 
ed tracheal branching during Drosophila development [1]. Multiple 
orthologues of spry have been reported in different organisms, 
including mouse, human, chicken, Xenopus and zebrafish [2-7]. 
The function of spry in the Lepidoptera, however, is not known. 
Mammalian genomes contain four Spry genes {Spry 1—4) encoding 
proteins (32-34 kDa) smaller than Drosophila melanogaster Spry 
(63 kDa) [8]. Drosophila Spry and vertebrate Spry proteins have a 
highly conserved C-terminal cysteine-rich region responsible for 
the membrane localization of Spry through palmitoylation [9] . A 
short region in the N terminus contains a conserved tyrosine 
residue, which mediates the interaction with its signaling molecules 
that contain Src-homology-2 domains [10-15]. 

Spry proteins are a major class of ligand-inducible inhibitors of 
RTK-dependent signaling pathways [16-17]. RTKs control a 
wide variety of processes, including proliferation, differentiation, 
migration and survival, in multicellular organisms [18-19]. In the 
RTKs- mitogen-activated protein kinase (MAPK) signaling 
pathway, the activated MAPKs phosphorylate and activate 
numerous target proteins, including transcription factors that 
regulate the expression of different genes [8,20-22]. The results of 



earlier genetic experiments indicated that the inhibitory activity of 
Spry is upstream of the extracellular signal-regulated kinase (ERK) 
and downstream of the RTK [8]. Later studies suggested the 
precise point at which Spry intercepts RTK signaling varies 
depending on the biological context. Studies with Drosophila 
indicated that during eye development. Spry inhibits signaling 
downstream of the epidermal growth factor receptor (EGFR) and 
upstream of rous sarcoma (Ras) [1] but functions at the level of 
rapidly accelerated fibrosarcoma (Raf) during wing and ovary 
development [23]. 

RTKs-mediated signaling events must be regulated precisely 
both spatially and temporally to achieve refinement of an 
appropriate biological outcome [24-27]. A salient feature of the 
RTK signaling pathway is the transcriptional induction of negative 
regulators by the pathways that are eventually inhibited, thereby 
providing an effective mechanism for the coordination of signaling 
input with the physiological response [28-34] . One such negative 
regulator is Spry, a multifaceted negative-feedback repressor of 
RTK signaling in vertebrates and invertebrates [35-36]. Activa- 
tion of RTK leads to the phosphohpid-dependent translocation of 
Spry to the plasma membrane, where it is tyrosine phosphorylated 
by an Src-like kinase activity [35,37]. Spry terminates this pathway 
by inhibiting the activation of Ras. And the study of Ras is well 
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done in silkworm[38-42]. Unphosphorylated Spry might also 
block the Ras-ERK pathway by inhibiting Raf 1 activation through 
an independent mechanism [12]. At the transcription level, 
activation of RTK leads also to the expression of Spry [43] . 

The available experimental evidence points to several mecha- 
nisms, all of which involve the interaction of Spry with essential 
elements of the RTK-Ras-ERK/MAPK cascade. The MAPK 
cascade, one of the most ancient and evolutionarily conserved 
signaling pathways, is important for many processes in immune 
responses [44] . MAPKs often have crucial roles in virus infection. 
Viral infection disorders the normal host cellular pathways, some 
of which reflect the organism's response to infection, whereas 
others are due to modification of the cellular environment by the 
virus [45-48]. For instance, modulation of MAPK pathways is 
essential for infection and replication of hepatitis B, Epstein-Barr, 
adenovirus and vaccinia viruses [49-52]. Further, it has been 
shown that activation oi B. mori MAPKs BmERK and BmJNK are 
required for B. mori nucleopolyhedrovirus (BmNPV) infection in 
BmN cells [53]. 

We cloned and identified a homologue of D. melanogaster Spry 
from the B. mori genome, and named it BmSpry. Using RNAi and 
over-expression in vivo and in vitro, we proved that BmSpry has a 
function in antiviral defense through regulation of the activation of 
ERK. This is the first report that Spry protein is involved in the 
antivirus response in the Lepidoptera. 

Materials and Methods 

Silkworm strain, cell lines and viruses 

B. mori T)Z, SN and Nm DZ lines were from the Gene Resource 
Library of Domestic;ated Silkworm (Southwest University, China). 
The BmE cell line[54] was cultured at 27°C in GRACE medium 
supplemented with 10% (v/v) fetal bovine serum (FBS). The 
BmN4-SIDl cell line was cultured at 27°C in IPL-41 medium 
supplemented with 10% (v/v) FBS [55]. BmNPV (Guangdong 
strain, China) and BmNPV-GFP were used in this study. Viruses 
were propagated in BmE cells and silkworm larvae, and BV titers 
were determined by plaque assay [56] . 

The mortality of DZ SN and Nm DZ hnes after oral inoculation 
with wild type BmNPV per os of the newly exuviated 2nd or 4th 
instar larvae were measured as described [57-58]. 

cDNA cloning, RT-PCR and qPCR analysis of BmSpry 

Total RNA was prepared using the Total RNA Kit II (Promega) 
[59] . Total RNA (2 (ig) was reverse transcribed as described [59] . 
For cDNA cloning, the GC buffer set (TaKaRa) and LA taq 
(TaKaRa) were used for PGR amplification with G+C-rich genes 
[60-61]. The primers used were: BmSpryF 5'-CGGTCGT- 
TTGGTTGGAGG-3' and BmSpryR 5'-ATTTGCGTAGAAT- 
CCAAAATACTAAC-3'. The nucleotide sequence reported in 
this paper has been submitted to the DDBJ/EMBL/GenBank 
Database under accession number KJ736835 [BmSpry). The 
qPCR/RT-PCR primers were: BmSpryF 5'-GGAGTGAATG- 
CGCACGAGTT-3' and BmSpryR 5'-CGCTGAAGTTGTT- 
CGTTGATTTC-3'. The housekeeping gene BGIBMGA003186- 
TA was used as an internal control to standardize the variance 
among the different templates as described [57-58]. QPCR was 
done with the ABI StepOnePlus Real-Time PGR System and ABI 
7500 Fast Real-Time PGR System (Apphed Biosystems). 

Knockdown of BmSpry in BmE cells, BmN4-SID1 cells and 

in individuals 

The dsRNAs for BmSpry and DsRed were generated by using a 
RiboMAX Large Scale RNA Production System-T7 kit (Promega) 



[62]. The primers were: T7-BmSpryF 5'-TAATACGACTCAC- 
TATAGGGGTGAAGGCGGTGTTGTAGGA-3', T7-BmSpryR 
(5'-TAATACGACTCACTATAGGGAGCGATCACTGAGTG- 
GACGC-3', T7-DsRedF 5'-TAATACGACTCACTATAGGG- 
GTACGGCTCCAAGGTGTACG-3', T7-DsRedR 5'-TAATA- 
GGAGTGACTATAGGGGGTGTAGTCGTCGTTGTGGG-3'. 

BmE cells (5x10^) were transfected with 4 \ig of dsRNA using 
the TransMessenger'^ Transfection Reagent (Qiagen) [63]. At 3 
days post-transfection, cells were challenged with BmNPV-GFP at 
an MOI of 1. At 3 days after infection, cells were processed for 
immunofluorescence and harvested for qPCR [64]. 

BmN4-SIDl cells (5x10"^) were used for RNAi by adding 
500 ng of dsRNA directly into 1 ml of medium supplemented with 
10% FBS [55]. After 5 days of RNAi, the rest of the process was as 
described above for the BmE cells. 

Total DNA was obtained from the cells with the Universal Genomic 
DNA Extraction Kit Ver 3.0 (TaKaRa) as described [65]. The qPGR 
primer GP64 was used for analysis of BmNPV replication with 
primers: GP64F 5'-CCATCGTGGAGACGGACTA-3', GP64R 5'- 
GTCGGACTGCTGCCTGA-3', BmGAPDHF 5'-CATTGGGG- 
GTGGGTGTTGGTAAT-3', BmGAPDHR 5'-GGTGGCTGC- 
TTGACCTTTTGC-3'. The housekeeping gene BmGAPDH was 
used as an internal control to standardize the variation among the 
dilferent templates. 

Nm DZ newly exuviated 5th instar larvae were injected with 
30 (ig of dsRNA [66] . Three days after RNAi, the larvae were 
injected with 2 |a.l of virtis (10^ pfu/ml) by stab inoculation as 
described [67]. Total DNA was obtained at 3 days after infection. 
Each sample was extracted from five treated larvae and used for 
qPCR analysis of BmNPV replication as described [58]. 

Over-expression of BmSpry in BmE cells 

The BmSpiy gene was amplified from the Nm DZ cDNA. The 
termination signal SV40 was cloned from the piggyBac [3 xp3 
EGFP afm] vector. The BmActin 4 (A4) promoter was used for the 
vector. The three elements were added to the empty vector 1 180. 
BmE cells (5x10^) were transfected with 0.8 ng of [A4-BmSpry- 
SV40] vector using the X-tremeGENE HP DNA Transfection 
Reagent (Roche) [68]. At 2 days after transfection, cells were 
challenged with BmNPV-GFP at an MOI of 1. At 3 days after 
infection, cells were processed for immunofluorescence and 
harvested for qPCR as described [64]. 

Western blotting 

BmE cells (5x10 ') were transfected with 0.8 |ig of vector. At 2 
days after transfection, cells were challenged with BmNPV-GFP at 
an MOI of 10. After infection for 24 h the cells were harvested for 
western blotting. RIPA lysis buffer was used to extract the total 
proteins within PhosSTOP (Roche) and PMSF (Roche) and 10 |J,g 
of total protein was used for western blotting analysis as described 
[53]. The antibody of anti-phospho-ERK (Cell Signaling Tech- 
nology) and anti-GAPDH (Sigma Aldrich) were used in this 
experiment. 

BmN4-SIDl cells (5x10^) were added to 500 ng of dsRNA in 
1 ml of meditim. After 5 days of RNAi, cells were challenged with 
BmNPV-GFP at an MOI of 10. After infection for 24 h, tiie cells 
were harvested for western blotting as described [53]. 

Results 

Characteristics of the BmSpry sequence 

To analyze the features of the BmSpiy sequence, we cloned this 
gene using LA taq and GC buffer. The BmSpry cDNA sequence 
contains 1398 nucleotides (Fig. 1) and the open reading frame. 
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which consists of 642 nucleotides, is G+C-rich (71.2%). The 
sequence encodes a propeptide consisting of 213 amino acid 
residues, which has a calculated molecular mass of 22,717 Da. 
The predicted isoelectric point (pi) of the mature peptide is 8.90. 

BmSpry, similar to D. melanogaster Spry, is a single-copy gene. The 
predicted protein sequence contains a Sprouty domain in its C 
terminus, which is conserved from invertebrates to vertebrates. 
This cysteine-rich domain consists of 104 amino acids (positions 
100-203). Although the N terminus was more variable than the C 
terminus, the motif at amino acid positions 45-5 1 contained an 
invariant tyrosine residue (^^47) conserved between invertebrates 
and vertebrates (Fig. 1). Many of the inhibitory functions of Spry 
proteins are dependent on this tyrosine residue [43 ,69]. 

Expression pattern of BmSpry 

The expression profile of BmSpry during the developmental 
stages from embryo to adult was investigated by RT-PCR using 
cDNAs from different tissues as templates. BmSpry was transcribed 
in all of the tissues examined (Fig. 2A) and the expression level of 
BmSpry was similar in all of these tissues, except the integument 
(Fig. 2A, lane 2). BmSpry was transcribed in all developmental 
stages from embryo to adult, but the expression levels were 
markedly higher in the egg, pupa, 4th and 5th instar compared to 
the other stages of development (Fig. 2B). 

Next, we used time-course analysis in the BmE cells by 
quantitative PGR (qPCR) to analyze the induced expression level 
of BmSpry after BmNPV infection (Fig. 2C). We chose five post- 
infection time points as the experimental group and time zero as a 
negative control. The results indicated the BmSpry expression level 
was down-regulated gradually from 2-24 h and was down- 
regulated markedly by BmNPV infection at 1 2 h and 24 h. These 
results implied BmSpry might have a function in the process of 
BmNPV infection. 

Restricting BmNPV infection in cultured cells 

To determine whether BmSpry was involved in antiviral defense 
against BmNPV, we generated the double-stranded RNA (dsRNA) 
oi BmSpry and repressed BmSpry in BmE cells via RNA interference 
(RNAi). The results showed silencing for BmSpry was efficient and 
the expression level of BmSpry was markedly reduced (Fig. 3A). 
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Figure 2. Expression profile of BmSpry in different tissues, 
developmental stages and inducible expression by BmNPV. (A) 

Expression of BmSpry in multiple tissues of 3^''-day 5* instar larvae. 
Lanes: 1, head; 2, integument; 3, hemocyte; 4, malpigiiian tubule; 5, 
midgut; 6, fat body; 7, silk gland; 8, testis; and 9, ovary. (B) Expression of 
BmSpry at developmental stages from embryo to adult. Lanes: 1 , 2 days 
after egg laying (AEL); 2, 4 days AEL; 3, 6 days AEL; 4, 8 days AEL; 5, 
newly hatched larvae; 6, molting larvae (ML) of the 1st instar; 7, newly 
exuviated larvae (NEL) of the 2nd instar; 8, ML of the 2nd instar; 9, NEL 
of the 3rd instar; 10, ML of the 3rd instar;1 1, NEL of the 4th instar; 12, ML 
of the 4th instar; 13, NEL of the 5th instar; 14, pupation in (PI) 2 days; 15, 
PI 4 days; 16, PI 6 days; 17, PI 8 days; and 18, adult moth. (C) Analysis of 
BmSpry gene expression in BmNPV-infected BmE ceils at MOI of 10. 
Total RNA was extracted from the BmE cells at the indicated time points 
post-infection. Non-infected cells were used as the control and qPCR of 
BmSpry was done. A representative of triplicate experiments is shown. 
Data are given as mean ±SD (n = 3). Statistically significant differences: 
** P<0.01. 

doi:1 0.1 371/journal.pone.0099200.g002 

Next, we challenged RNAi-treated cells with BmNPV-green 
fluorescent protein (GFP) and then investigated the BmE cells by 
fluorescence microscopy at 3 days post-infection. The result 
showed a significant increase of infected cells upon silencing 
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TACGAGTACCTTGTAAATAT AGTTTTATAGTGTGAATGT A TGTGAACGTGTCGCCCGAAC GtCGAAAGTCTCAAATTTTC TTGAGACCGTAACGTGTCAT 
AAAAAAAACGTTAACGTTAG TATTTTGGATTCTACGCAAA TCCAGAGAAGATTAAGAAAT GAGATATAAAGCATACTTAT TAAGATGTAGATAATTTACC 
ATTTATTTATTAACTGAATT TGTATTTATTTGTAAATAAA ACGACAACATTATATTATAA ATACGAATGTTTACTACATA aaaaaaaaaaaaaaaaaa 



Figure 1. Characteristics of the Fm5p/y sequence. The nucleotide sequence is 1398 bp long and the predicted protein contains 213 amino 
acids. Numbering of the DNA sequence in the left-hand margin begins with the A of the presumptive initiating AUG codon. The 213 amino acid 
sequence is presented, terminated by an asterisk (*). The N terminus has a short conserved motif (underlined) and a conserved residue Y (indicated by 
a circle). A Cys-rich region (21 residues) is marked C rich at the right. 
doi:1 0.1 371/journal.pone.0099200.g001 
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Figure 3. 5m5p/y inhibited BmNPV replication in cultured cells. (A) BmE cells pretreated with dsRNA as indicated, the dsRNA of dsRed was 
used as a negative control. At 3 days post transfection, total RNA was extracted and qPCR was used to analyze the BmSpry expression level. (B) BmE 
cells treated with dsRNA against the indicated genes were infected with BmNPV-GFP at IVIOI of 1 for 3 days and processed for immunofluorescence. 
(C) BmE cells treated with the indicated dsRNA were infected with BmNPV-GFP at IVIOI of 1 and infection total genomes were extracted for qPCR at 3 
days post. (D) BmE cells were used for transient transfection and the empty vector 1180 was used as a negative control. BmSpry-oe was an over- 
expression vector of BmSpry. At 3 days post-transfection, total RNA was extracted for qPCR. (E) BmE cells were subjected to transient transfection 
with BmSpry expression vectors and the empty 1 180 vector, as indicated. At 2 days post-transfection, the cells were infected with BmNPV-GFP at MOI 
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of 1 for 3 days and then processed for immunofluorescence. (F) BmE cells transfected with the indicated vectors were infected with BmNPV at MOI of 
1 and total genomes were extracted for qPCR at 3 days post-infection. A representative of triplicate experiments is shown. Data are given as mean 
±SD (n = 3). Statistically significant differences: ** P<0.01. 
doi:l 0.1 371/journal.pone.0099200.g003 



BmSpry but not the Discosoma red fluorescent protein (dsRed) 
(Fig. 3B). We extracted the total genomes of the infected cells for 
qPCR analysis and the results indicated there was a significant 
increase in the amount of viruses of about twofold compared to the 
BmE cells when depleting BmSpry by RNAi (Fig. 3C). 

We used BmN4-SID 1 cells to repeat the RNAi experiment for 
further verification of the RNAi results and obtained the same 
results as those found for the BmE cells. Firstiy, \vc confirmed the 
BmN4-SID 1 cells reduced the BmSpry expression level significantly 
without the use of a transfection reagent (Fig. SI A). Then, we 
challenged RNAi-treated BmN4-SIDl cells with BmNPV and 
analyzed the amount of viruses, as we did with BmE cells. The 
results showed the amount of viruses in the BmSpry-depleted cells 
was more than twofold greater compared to the control cells (Fig. 
SIB). 

We up-regulated the expression level of BmSpiy in the BmE cells 
by transient transfi'ction and ()\'er-expression of BmSpiy was 
verified by qPCR (Fig. 3D). At 48 h post transfection, we 
challenged over-expression-treated cells with BmNPV-GFP and 
then analyzed the infection by fluorescence microscopy. The 
percentage of infected cells was notably reduced in the cells with 
over-expression ot BmSpry (Fig. 3E). Furthermore, qPCR analysis 
confirmed over-expression ot BmSpry inhibited BmNPV replication 




Figure 4. BmSpry regulation of ERK activation inhibited 
BmNPV replication. (A) Firstly, BmE cells were transfected. The 
empty vector 1 180 was used as a negative control and the BmSpry-oe 
vector was used as the experimental group with over-expression of 
BmSpry. At 2 days post-transfection, BmE cells were infected witli 
BmNPV at MOI of 10. Post infection for 24 h, total proteins were 
extracted for western blotting using anti-phosplno-ERK antibody. Total 
protein levels were measured with anti-GAPDH antibody. (B) BmN4- 
SID1 cells were used. After RNAi of 5 days, the BmN4-SID1 cells were 
infected with BmNPV at MOI of 10. Post infection for 24 h, total proteins 
were extracted for western blotting. A representative of triplicate 
experiments is shown. 
doi:1 0.1 371/journal.pone.0099200.g004 



and reduced the amount of viruses to half (Fig.SF). Taken 
together, our data suggest that BmSpry has an important role in 
BmNPV infection in vitro. 

BmSpry inhibits BmNPV replication through regulation of 
ERK activation 

Early work showed activation of ERK and JNK is required for 
BmNPV infection in cultured cells [5 3]. In this study, we used over- 
expression and RNAi of BmSpry in BmNPV-infected cells to 
investigate whether the antiviral defense function of BmSpry acts 
through inhibition of ERK activation. ERK had a basic level of 
phosphorylation in uninfected cells. When the cells had been 
infected with BmNPV for 24 h, however, the activation of ERK 
was enhanced significandy. Interestingly, the activation state of 
ERK induced by BmNPV infection can be inhibited markedly in 
BmE cells by over-expression of BmSpry (Fig. 4A). Furtlu'rmore, 
down-regulation of BmSpiy in BmN4-SID 1 cells can enhance the 
activation state of ERK significandy by BmNPV infection (Fig. 4B). 
These results suggested regulation of the ERK pathway by the 
BmSpry protein inhibited replication of the virus. 

Antiviral defense in vivo 

We showed BmSpry was antiviral in vitro and asked whether it 
had similar biological function in vivo. Using RNAi in vivo, we 
investigated whether loss of BmSpry had an effect on BmNPV 
replication. Knock-down BmSpry in silkworm larvae was generated 
by injection of 30 )j.g of dsRNA against BmSpry. Efficient silencing 
\vas confirmed by qPCR; the expression level of BmSpry was 
reduced by 40% (Fig.5A). Subsequentiy, we challenged RNAi- 
treated larvae with BmNPV by stab inoculation and analyzed viral 
DNA replication with qPCR. There was a significant increase, 
about twofold, of viral DNA replication upon silencing oi BmSpry 
but not the dsRed (Fig. .5B). 

Next, we used silkworm lines DZ SN and Nm DZ to investigate 
the mortality following infection with BmNPV per os. The Nm DZ 
line was the wild type and the DZ SN line had the same genetic 
background apart from a transposon insertion 14.6 kb upstream of 
BmSpiy. Interestingly, the expression level BmSpry was markedly 
lower in DZ SN compared to Nm DZ (Fig. 5C). Newly exuviated 
4th-instar larvae of the two lines were infected with BmNPV per os 
using a dose of 10* or lO*" occlusion bodies (OB)/larva. The results 
showed mortality was significandy higher in the DZ SN line 
compared to the Nm DZ line for each dose of OB. At the lower 
dose, mortality was nearly 10% in the Nm DZ line and reached 
60% in the DZ SN line. At the higher dose, mortality was 50% in 
tiie Nm DZ line and 90% in tiie DZ SN line (Fig. 5D). 

The same results for mortality were found in newly exuviated 
2nd-instar larvae following infection with BmNPV per os at a dose 
of 5 X 10^ OB/larva. The Nm DZ line had a mortality of —5% and 
the DZ SN fine had a mortality >20% (Fig. 5E). Overall, these 
results showed the DZ SN line was significandy susceptible to 
BmNPV, suggesting BmSpry has an important antiviral role in vivo. 

Discussion 

The essential roles for two MAPKs, ERK and JNK, in the 
efficient infection of silkworm cells was firmly established [53]. It 
was unclear, however, how silkworm larvae resist virus infection 
using these signaling pathways and how viruses overcome the host 
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Figure 5. BmSpry^Nas essential for antiviral defense in vivo. (A) The silkworm Nm DZ line was used for the RNAi experiment. Injection of dsRNA 
is indicated, the dsRNA of dsRed was used as a negative control with dsRNA of BmSpry as the experimental group. At 3 days post-injection of dsRNA, 
total RNA was extracted for qPCR to measure the expression level of BmSpry. Data are given as mean ±SD (n = 3). (B) Nm DZ silkworms treated with 
dsRNA against the indicated genes were infected with BmNPV (10^ pfu/mL) by stab inoculation for 3 days and processed for qPCR. Data are given as 
mean ±SD (n = 3). (C) Analysis of the expression level of BmSpry in Nm DZ and the mutant DZ SN. Total RNA extracted from newly exuviated Nm DZ 
and DZ SN 5th instar larvae was used for qPCR. Data are given as mean ±SD (n = 3). (D) Analysis of mortality of DZ SN and Nm DZ after oral 
inoculation with BmNPV. Nm DZ and DZ SN newly exuviated 4th instar larvae were used for this experiment. DZ SN and Nm DZ silkworms were 
infected with BmNPV per os using a dose of lO'' or 10* OB/larva and mortality was monitored until the adult stage (22 days). (E) Newly exuviated 2nd- 
instar larvae were used to investigate mortality after inoculation per os with 5x10^ OB/larva. Mortality was monitored until the pupa stage (16 days). 
A representative of triplicate experiments is shown. Statistically significant differences: ** P<0.01. 
doi:1 0.1 371/Journal.pone.0099200.g005 



resistance. In this study, we cloned a gene and named it BmSpry. 
Expression pattern analysis showed it was transcribed in all larval 
tissues and all developmental stages from egg to adult moth. We 
showed also BmSpry had the ability to normalize the activation of 
ERK induced by BmNPV infection and the virus had the ability to 
down-regulate BmSpry, rendering the MAPK signaling pathway 
out of control. This is the first report of this important function of 
BmSpry in antiviral defense. 

SpTj), a general inhibitor of RTKs [16-17], has many copies in 
some organisms; e.g. four copies in the human and mouse 
genomes [8]. There is only one copy in the silkworm, however, 
which is likely why BmSpry was expressed in all of the larval tissues 
and at all developmental stages [70-72]. The expression level of 
BmSpyy was very similar among all tissues, suggesting BmSpry has 
the same function in all silkworm tissues. At different develop- 
mental stages, however, there was a great diversity of expression 
level, with especially high levels in the egg and pupa. This result 
implies BmSpry has a crucial role in development of the individual. 
The spry gene must regulate RTK-mediated signaling pathways 
precisely to ensure the appropriate biological outcome. 



It is proven that spry is a general inhibitor of RTKs [16-1 7] and 
each RTK signaling pathway has a salient feature requiring exact 
feedback regulation systems [28—34]. In a normal individual, 
activation of ERK is foUo^vcd by translocation to the nucleus and 
induction of spry transcription [43]. Feedback of the newly 
synthetized spry protein regulates the MAPK signaling pathway 
to the normal level. In the silkworm, stable expression of BmSpry 
ensures the appropriate biological outcome. When the virus 
disequilibrates the RTK-AIAPK signaling pathway and enhances 
the level of ERK activation, BmSpry acts as an RTK signaling 
pathway inhibitor restoring the balance and restricts the use of 
RTKs by BmNPV for replication. 

Baculovirus-host interactions are a hot topic for research but 
progress is slow. Since it was found baculovirus use the MAPK 
signaling pathways for replication [53], the search for the receptor 
of the upstream of MAPK has never been abandoned. In this 
study, we showed BmSpry is involved in the MAPK signaling 
pathway, indicating the receptor of the upstream of MAPK 
belongs to the RTK family. This finding allows the search for the 
receptor to be more explicit and specific. Further, we suggest the 
baculovirus uses the host growth factors to activate the RTKs 
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because the BmNPV genome contains only one of the growth 
factors, vFGF [73], and is not necessary for BmNPV to activate 
MAPK [53]. In the BmNPV-infected ceUs, the virus could activate 
the ERK and down-regulate BmSpry simultaneously. This infection 
mechanism allows the virus unlimited use of the cell resources for 
maximum replication. 

In conclusion, our results show that BmSpry is involved in 
BmNPV infection. BmSpry is able to act against baculoviruses, 
modulating the MAPK signaling pathways and normalizing the 
superactivation of ERK by virus infection. The results of our study 
contribute to elucidation of the mechanism that allows baculovirus 
to modulate the MAPKs and confirm RTKs are upstream of 
MAPKs during baculovirus infection. 

Supporting Information 

Figure SI BmSpry inhibited BmNPV replication in BmN4-SIDl 

cells. (A) The BmN4-SIDl cells were used for the RNAi 
experiment and the dsRNA of dsRed was used as a negative 

References 

1. Hacohen N, Kramer S, Sutherland D, Hiromi Y, Krasnow MA (1998) sprouty 
Encodes a Novel Antagonist of FGF Signaling that Patterns Apical Branching of 
the Drosophila Airways. Cell 92: 253-263. 

2. Chambers D, Medhurst A, Walsh F, Price J, Mason I (2000) Differential Display 
of Genes Expressed at the Midbrain-Hindbrain Junction Identifies sprouty2: An 
FGF8-Induciblc Member of a Family of Intracellular FGF Antagonists. 
Molecular and Cellular Ncuroscicncc 15: 22-35. 

3. dc Maximy AA, Nakatakc Y, Moncada S, Itoh N, Thicry JP, et al. (1999) 
Cloning and expression pattern of a mouse homologue of Drosophila sprouty in the 
mouse embr\'o. Mechanisms of development 81: 213—216. 

4. Furthauer M, Reifers F, Brand M, Thisse B, 'Fhisse C (2001) sprouty4 acts in vivo 
as a feedback-induced antagonist of F(iF signaling in zebralish. Development 
128: 2175-2186. 

5. Minowada G, Jarvis LA, Chi CL, Neubuscr A, Sun X, ct al. (1999) Vertebrate 
Sprouty genes are induced by FGF signaling and can cause chondrodysplasia 
when overexpressed. Development 126: 4465^475. 

6. Nutt SL, Dingwell KS, Holt CE, Amaya E (2001) Xenopus Sprouty2 inhibits 
FGF-mediaicd gastrulation movements but does not affect mesoderm induction 
and patterning. (Jenes & development 13: I 152 1 166. 

7. TelftJD, Lee M, Smith S, Leimvand M. /IkkiJ. ci al. ;1999) Conscrxrd function 

mSpry-2, a murine homolog of Dro.^opliila sprouty, which negatively modulates 
respiratory organogenesis. Current biolog\' 9: 219-222. 

8. Kim HJ, Bar-Sagi D (2004) Modulation of signalling by Sprouty: a developing 
story. Nature Reviews Molecular Cell Biology 5: 441-450. 

9. Sasaki A, Taketomi T, Kato R, Saeki K, Nonami A, et al. (2003) Mammahan 
Sprouty4 suppresses Ras-independent ERK activation by binding to Rafl. 
Nature ceU biology 5: 427^32. 

10. Fong CW, Lcong HF. Wong ESM, Lim J, Yusoff P, et al. (2003) Tyrosine 
phosphorvlalion of Sproiity2 enhances its interaction with c-Cbl and is crucial 
for its function. Journal of Biological Chemistry 278: 33456-33464. 

11. Hall AB, Jura N, DaSilvaJ, Jang YJ, Gong D, et al. (2003) hSpr\'-2 is targeted to 
the ubiquitin-dependent proteasome pathway by c-Cbl. Current biology 13: 
308-314. 

12. Hanafusa H, Torii S, Yasunaga T, Nishida E (2002) Sproutyl and Sprouty2 
provide a control mechanism for the Ras/MAPK signalling pathway. Nature 
ceU biology 4: 850-858. 

13. Mason JM, Morrison DJ, Bassit B, Dimri M, Band H, et al. (2004) Tyrosine 
phosphorylation of Sprouty proteins regulates their ability to inhibit growth 
factor signaling: a dual feedback loop. Molecular biology of the cell 15: 2176— 
2188. 

14. Rubin C, Litvak V, Mcdvedovsky H, Zwang Y, Lev S, et al. (2003) Sprouty fme- 
tuncs EGF signaling through interlinked positive and negative feedback loops. 
Current biolog>' 13: 297-307. 

15. Smart J, Oppermann H, Czernilofsky A, Purchio A. Erikson R, ct al. (1981) 
Characterization of sites for tyrosine phosphorylation in the translbrming protein 
of Rous sarcoma virus (pp60v-src) and its normal cellular homologue (pp60c-src). 
Proceedings of the National Academy of Sciences 78: 6013-6017. 

16. Gross I, Bassit B, Benezra M, Licht JD (2001) Mammalian sprouty proteins 
inhibit cell growth and differentiation by preventing ras activation. Journal of 
Biological Chemistry 276: 46460-46468. 

17. Yusoff P, Lao D-H, Ong SH, Wong ESM, Lim J, et al. (2002) Sprouty2 inhibits 
the Ras/MAP kinase pathway by inhibiting the activation of Raf Journal of 
Biological Chemistry' 277: 3195 3201. 

18. Blume-Jensen P, Hunter T (2001) Oncogenic kinase signaUing. Nature 411: 
355-365. 



control. After 5 days of RNAi, total RNA was extracted for qPCR. 
(B) BmN4-SIDl cells treated with the indicated dsRNA were 
infected at MOI of 1 and the total genomes were extracted for 
qPCR at 3 days post- infection. A representative of triplicate 
experiments is shown. Data are given as mean ±SD (n = 3). 
Statistically significant differences: ** P<0.01. 
(TIF) 

Acknowledgments 

We thank Dr. Xiaofeng Wu for gifting the BmNPV-GFP virus. 

Author Contributions 

Conccix'Cfl and designed the experimenls: SJ TC QX. Performed the 
experiments: SJ LJ QY YX. Analyzed the data: SJ TC PL. Contributed 
reagents/materials/analysis tools: QY YX TK QX. Wrote the paper: SJ 
TC QX. 



19. Ullrich A, Schlessingcr J (1990) Signal transduction by receptors with tyrosine 
kinase activity. Cell 61: 203-212. 

20. Impagnatiello M-A, Weitzer S, Gannon G, Compagni A, Gotten M, et al. (200 1) 
Mammalian sprouty- 1 and-2 are membrane-anchored phosphoprotein inhibi- 
tors of growth factor signaling in endothelial cells. The Journal of cell biology 
152: 1087-1098. 

21. Ozaki K, Kadomoto R, Asato K, Tanimura S, Itoh et al. (2001) ERK 
pathway positively regulates the expression of Sprouty genes. Biochemical and 
biophvsical research communications 285: 1084—1088. 

22. Sasaki A, Taketomi T, Wakioka T, Kato R, Yoshimura A (2001) Identification 
of a dominant negative mutant of Sprouty that potentiates fibroblast growth 
factor-but not epidermal growth factor-induced ERK activation. Journal of 
Biological Chemistry 276: 36804-36808. 

23. Reich A, Sapir A, Shilo B (1999) Sprouty is a general inhibitor of receptor 
tyrosine kinase signaling. Development 126: 4139—4147. 

24. Lax I, Wong A, Lamothe B, Lee A, Erost A, et al. (2002) The docking protein 
FRS2c( controls a MAP kinase-mediated negative feedback mechanism for 
signaling by E(iE receptors. Molecular cell 10: 709—719. 

25. Ostman A, Bohmer E-D (2001) Regulation of receptor tyrosine kinase signaling 
by protein tyrosine phosphatases. Trends in cell biology 11: 258—266. 

26. Reynolds AR, Tischer C, Verveer PJ, Rocks O, Bastiaens PI (2003) EGFR 
activation coupled to inhibition of tyrosine phosphatases causes lateral s^al 
propagation. Nature cell biology 5: 447-453. 

27. Tonks NK (2006) Protein tyrosine phosphatases: from genes, to function, to 
disease. Nature Reviews Molecular Cell Biology 7: 833—846. 

28. Amit I, Citri A, Shay T, Lu Y, Katz M, et al. (2007) A module of negative 
feedback regulators defines growth factor signaling. Nature genetics 39: 503- 
512. 

29. Buday L, Warne PH, DownwardJ (1995) Downregulation of the Ras activation 
pathway by MAP kinase phosphorylation of Sos. Oncogene 11: 1327-1331. 

30. Gual P, Giordano S, Anguissola S, Parker P, ComogHo P (2001) Gabl 
phosphorylation: a novel mechanism for negative regulation of HGF receptor 
signaling. Oncogene 20: 156—166. 

31. Mothe I, Van Obberghen E (1996) Phosphorylation of insulin receptor 
substrate-1 on multiple serine residues, 612, 632, 662. and 731, modulates 
insulin action. Journal of Biologic'al (.:hrmislry 271: 1 1222-11227. 

32. Shilo B-Z (2005) Regulating ihr dynamics of E(T'' receptor signaling in space 
and time. Development 132: 4017 4027. 

33. Sorkin A, Goh LK (2009) Endocytosis and intracellular trafficking of ErbBs. 
Experimental cell research 315: 683-696. 

34. Ueki K, Matsuda S, Tobe K, Gotoh Y, Tamemoto H, et al. (1994) Feedback 
regulation of mitogen-activated protein kinase kinase kinase activity of c-Raf-1 
by insulin and phorbol ester stimulation. Joumail of Biological Chemistry 269: 
15756-15761. 

35. Casci 1\ Vinos J, Freeman M (1999) Sprouty, an intracellular inhibitor of Ras 
signahng. Cell 96: 655-665. 

36. Kramer S, Okabe M, Hacohen N, Krasnow MA, Hiromi Y (1999) Sprouty: a 
common antagonist of FCiF and EGF signaling pathways in Drosophila. 
Development 126: 2515-2525. 

37. Lim J, Wong ESM, Ong SH, Yusoff P, Low BC, et al. (2000) Sprouty Proteins 
Are Targeted to Membrane Ruffles upon Growth Factor Receptor Tyrosine 
Kinase Activation IDENTIFICATION OF A NOVEL TRANSLOCATION 
DOMAIN. Journal of Biological Chemistry 275: 32837-32845. 

38. Ma L, Liu S, Shi M, Chen X-x, Li S (2013) Ras 1 CA-upregulated bcpi inhibits 
cathepsin activity to prevent tissue destruction of the Bombyx posterior silk 
gland. Journal of proteome research 12: 1924—1934. 



PLCS ONE I www.plosone.org 



7 



June 2014 I Volume 9 | Issue 6 | e99200 



BmSpry Resist BmNPV Reproduction 



39. Ma L, Ma Li X, Cheng L, Li K, ct al. (2014) Transcriptomic analysis of 
differentially expressed genes in the Ras 1 C A-overexpressed and wildtype 
posterior silk glands. BMCJ genomies 15: 182. 

40. Ma L, Xu H, Zhu J, Ma S, Liu Y, ct al. (201 1) RaslCA overexpression in the 
posterior silk gland improves silk yield. Cell research 21: 934—943. 

41. Moriya K, Tsubota T, Ishibashi N, Yafime A, Suzuki T, et al. (2010) Bombyx 
mori Ras proteins BmRasl, BmRas2 and BmRas3 are neither farnesylated nor 
palmitoylatcd but are gcranylgcranylated. Insect molecular biolo^g)' 19: 291-301 . 

42. Ogura T, Tan A, Tsubota T, Nakakura T, Shiotsuki T (2009) Identification and 
expression analysis of ras gene in silkworm, Bombyx mori. PloS one 4: e8030. 

43. Mason JM, Morrison DJ, .\lbert Basson M, Licht JD (2006) Sprouty proteins: 
multifaceted negati\'e-feedback regulators of receptor tyrosine kinase signaling. 
Trends in cell biology 16: 45-54. 

44. Dong C, Davis RJ, Flavell RA (2002) MAP kinases in the immune response. 
Aimual review of immunology 20: 55—72. 

45. Ludwig S, Planz O, Pleschka S, Wolff T (2003) Influenza-virus-induced signaling 
cascades: targets for antiviral therapy? Trends in molecular medicine 9: 46-52. 

46. Panteva M, Korkaya H, Jamecl S (2003) Hepatitis viruses and the MAPK 
pathway: is this a survival strate^gy? Virus research 92: 131-140. 

47. Tyler KL, Clarke P, DcBiasi RL, Kominsky D, Poggioli GJ (2001) Reoviruses 
and the host cell. TRENDS in Microbiology 9: 560-564. 

48. Xiao W, Yang Y, Weng Lin T, Yuan M, et al. (2009) The role of the PI3K- 
Akt signal transduction pathway in Autographa califomica multiple nucleopolyhedrovirus 
infection oi Spodoptera frugiperda cells. Virology 391: 83—89. 

49. de Magalhacs JC, Andrade AA, Silva PN, Sousa LP, Ropert C, et al. (2001) A 
Mitogenic Signal Triggered at an Early Stage of Vaccinia Virus Infection 
IMPUCATION OF MEK/ERK AND PROTEIN KINASE A IN VIRUS 
MULTIPLICATION. Journal of Biological Chemistry 276: 38353-38360. 

50. Gao X, Wang H, Sairenji T (2004) Inhibition of Epstein-Barr virus (EBV) 
reactivation by short interfering RNAs targeting p38 mitogen-activated protein 
kinase or c-myc in EBV-positive epithelial cells. Journal of virology 78: 1 1798— 
11806. 

51. Ikeda M, Kobayashi M (1999) Cell-Cycle Perturbation in Sf9 Cells Infected with 
Autographa ca^mico Nucleopolyhedrovirus. Virology 258: 176-188. 

52. Zaremba A, Schmuecker U, Esche H (201 1) Sprouty is a cytoplasmic target of 
adenoviral El A oncoproteins to regulate the receptor tyrosine kinase signalling 
pathway. Virology journal 8: 192. 

53. Katsuma S, Mita K, Shimada T (2007) ERK-and JNK-dependent signaling 
pathways contribute to Bombyx mori nucleopolyhedrovirus infection. Journal of 
virology' ai: 13700 -13709. 

54. Pan M-H, Xiao S-Q, C.:hcn M, Hong X-J, Lu C (2007) Establishment and 
characterization of two embryonic cell lines of Bombyx mori. In Vitro OUular & 
Developmental Biology -^Animal 43: 101-104. 

55. Mon H, Kobayashi I, Ohkubo S, Tomita S, Lee J, et al. (2012) Effective RNA 
interference in cultured silkworm cells mediated by overexpression of 
Caenorhabditis elegans SID-1. RNA biology 9: 40--1-6. 

56. Katsuma S, Noguchi Y, Zhou C, Kobayashi M, Maeda S (1999) Character- 
ization of the 25K FP gene of the baculovirus Bombyx mori nucleopolyhe- 
drovirus: implications for post-mortem host degradation. Journal of General 
Virology 80: 783-791. 

57. Jiang L, Cheng T, Zhao P, Yang Q, Wang G, et al. (2012) Resistance to 
BmNPV via Overexpression of an Exogenous Gene Controlled by an Inducible 



Promoter and Enhancer in Transgenic Silkworm, Bombyx mori. PloS one 7: 

e41838. 

58. Jiang L, Wang G, Cheng T, Yang Q, Jin S, et al. (2012) Resistance to Bombyx 
mori nucleopolyhedrovirus via overexpression of an endogenous antiviral gene 
in transgenic silkworms. Archives of virology 157: 1323-1328. 

59. Mortensen AS, Arukwe A (2006) Dimethyl sulfoxide is a potent modulator of 
estrogen receptor isofbrms and xenoestrogcn biomarker responses in primary 
culture of salmon hepatocytes. Ac^uatic toxicology 79: 99-103. 

60. jiano- X, On M, Dcnholm I, Fang J, Jiang \V, ct al. (2009) Mutation in 
a( ii\ Irhnlincslrrase I associared wiih iriazophos resistance in rice stem borer, 
Chilo sappressalis (Lepidoptcra: Pyralidac). Biochemical and biophysical research 
communications 378: 269-272. 

61. Soejima K, Mimura N, Hiroshima M, Maeda H, Hamamoto T, et al. (2001) A 
Novel Human Metalloprotease Synthesized in the Liver and Secreted into the 
Blood: Possibly, the von WUlebrand Factor — Cleaving Protease? Journal of 
biochemistry 130: 475-480. 

62. Payungporn S, Chutinimitkul S, Ghaisingh A, Damrongwantanapokin S, 
Buranathai G, et al. (2006) Single step multiplex real-time RT-PCR for H5N1 
influenza A virus detection. Journal of virological methods 131: 143-147. 

63. Krysan K, Dalwadi H, Sharma S, Pold M, Dubinett S (2004) Cyclooxygenase 2- 
dependent expression of survivin is critical for apoptosis resistance in non-smaU 
cell lung cancer. Cancer reseairch 64: 6359-6362. 

64. Nakamoto M, Moy RH, Xu J, Bambina S, Yasunaga A, et al. (2012) Virus 
Recognition by Toll-7 Activates Antiviral Autophagy in Drosophila. Immunity 36: 
658-667. 

65. Bao Y-Y, Tang X-D, Lv Z-Y, Wang X-Y, Tian G-H, et al. (2009) (;ene 
expression profiling of resistant and susceptible Bombyx mori strains reveals 
nucleopolyhedrovirus-associated variations in host gene transcript levels. 
Genomics 94: 138-145. 

66. Deng H, ZhangJ, Li Y, Zheng S, Liu L, et al. (2012) Homeodomain POU and 
Abd-A proteins regulate the transcription of pupal genes during metamorphosis 
of the silkworm, Bombyx mori. Proceedings of the National Academy of 
Sciences 109: 12598-12603. 

67. Zhang Y, Lv Z, Chen J, Chen Quan Y, et al. (2008) A novel method for 
isolation of membrane proteins: A baculovirus surface display system. 
Proteomics 8: 4178-4185. 

68. Ramirez CL, Certo MT, Mussohno C, Goodwin MJ, Cradick TJ, et al. (2012) 
Engineered zinc finger nickases induce homology-directed repair with reduced 
mutagenic effects. Nucleic acids research 40: 5560-5568. 

69. Ishida M, Ichihara M, Mii S, Jijiwa M, Asai N, et al. (2007) Sprouty2 regulates 
growth and differentiation of human neuroblastoma cells through RET tyrosine 
kinase. Cancer science 98: 815-821. 

70. Goldsmith MR, Shimada T, Abe H (2005) THE GENETICS AND 
GENOMICS OF THE SILKWORM, BOMBYX MORI*. Annu Rev Entomol 

50: 71-100. 

7 1 . Xia Quo Y, Zhang Z, Li D, Xuan Z, et al. (2009) Complete resequencing of 
40 genomes reveals domestication events and genes in silkworm (Bombyx). 
Science 326: 433-436. 

72. Xia Q, Zhou Z, Lu C, Cheng D, Dai F, et al. (2004) A draft sequence for the 

genome of the domesticated silkworm (Bombyx mori). Science 306: 1937-1940. 

73. Gomi S, Majima K, Alaeda S (1999) Sequence analysis of the genome of 
Bombyx mori nucleopolyhedrovirus. Journal of General Virology 80: 1323— 
1337. 



PLOS ONE I www.plosone.org 



8 



June 2014 I Volume 9 | Issue 6 | e99200 



